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Abstract 

We consider multiple lattices and functions defined on them. We in- 
troduce slow varying conditions for functions defined on the lattice and 
express the variation of a function in terms of an asymptotic expansion 
with respect to the slow varying lattices. 

We use these results to perform the multiple-scale reduction of the 
lattice potential Korteweg-de Vries equation. 

1 Introduction 

The reductive perturbation method (or multiple-scale analysis) [18] allows us 
to deduce a set of simplified equations starting from a basic model without 
loosing its main characteristic features. The method consists essentially in an 
asymptotic analysis of a perturbation series, based on the existence of different 
scales to cure secularity. 

The success of the method relies mainly on the nice property of the resulting 
reduced models, which are simple and often integrable. Simple here means ac- 
tually simpler than the starting equations and still providing useful information. 
Integrable means that they carry an infinite set of conserved quantities, have an 
infinite set of symmetries and of exact solutions. Finally, as emphasized in [2], 
the reductive perturbation approach preserve integrability. Consequently this 
approach can be used to obtain new integrable models from known ones. 

The situation is quite different in the case of differential equations on a lattice 
(for example, in the case of dynamical systems when one has a continuous time 
and discrete space variables) for which a reliable reductive perturbative method 
which would produce reduced discrete systems up to our knowledge does not 
exist. Leon and Manna [11] and later Levi and Heredero [12] proposed a set 
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of tools which allow to perform multiple-scale analysis for a discrete evolution 
equation. These tools rely on the definition of a large grid scale via the compar- 
ison of the magnitude of related difference operators and on the introduction 
of a slow varying condition for function defined on the lattice. Their results, 
however, are not very promising as the reduced models are neither simpler nor 
more integrable than the original one. Starting from an integrable model, like 
the Toda lattice [19], the Leon and Manna reduction technique produce a non- 
integrable differential difference equation of the discrete Nonlinear Schrodinger 
type [13,17]. Levi and Heredero [12] started from the integrable differential- 
difference Nonlinear Schrodinger equation and got a nonintegrable system of 
differential-difference equations of Kortewg-de Vries type. 

We consider here the case of completely discrete equations defined on a two 
dimensional orthogonal lattice. We follow the approach introduced by Levi and 
Heredero [12], extended to the case of multiple orthogonal lattices. We try to 
keep all passages consistent with the continuous limit, when the lattice spacings 
on the different grids go to zero. 

In Section |21 we introduce, following [12], the multiple lattices, the slow 
varying conditions and the asymptotic expansions of the functions' variations 
while in Section [31 we apply the resulting formulas to the case of the multiple- 
scale expansion of the lattice potential Korteweg - de Vries equation (lpKdV) 
[8,14], 

(p — q + U«,m+1 - Un+l,m)(P + Q~ "n+l.m+l + U n ,m) = p 2 - Q 2 ■ (1) 

At the end, in Section 0] we discuss the results obtained and present a list of 
open problems and remarks relevant also for the case of a differential-difference 
equation [12]. 

2 Multi-lattice structure and the variation of a 
function on them 

2.1 Rescaling on the lattice 

Given a lattice defined by a constant lattice spacing h, we will introduce an 
apriory infinite number of lattices defined by lattice spacings Hj, with j = 
1,2, ...,oo, where Hj are well defined functions of h, Hj = Hj(h). In Fig. 
^ we show an example of such a situation with j = 1,2. For convenience we 
will denote by rij the running index of the points separated by Hj and n those 
separated by h. Moreover, in correspondence with the lattice variables, we can 
introduce the real variables x = hn and Xj — Hjirij. 

A simple definition of Hj is obtained by introducing an integer number N 
and defining Hj — N^h. If TV is a large number than = e will be a small 
number. The variables x and Xj will go over to continuous variables when 
respectively h — » 0, n — ► oo and Hj — > 0, rij — > oo in such a way that their 
products x — nh and Xj — rijHj are finite. 
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Figure 1: Multiple lattices (n, m = jjin>2 = with N = 5; the corresponding 
continuous coordinate would be x — nh, x\ = niH\ and xi = 712-^2 where 
Hi = hN, H 2 = hN 2 . 

Let us assume that Xj = e?x. Then, if x ~ ^j, Xj ~ 1. So represents, as 
j increases an always larger portion of the x axis. This assumption, together 
with the choice e = -h, will reflect onto a relation between the lattice variables 
n and nj as 

Xj = HjUi = hN^n-i = e*x = — -hn => n., = 
j j j j jy 3 j 

Consequently we need to move N 2j points on the lattice of the discrete variable 
n to shift the lattice variable rij by 1 point. 

2.2 Slowly varying functions and their expansion 

Here we study the relation between functions acting on the different lattices 
defined in Section I^Tl 

Let us consider a function / defined on the points of a lattice of index n given 
in Fig. ^ i.e. f(n). We are interested in understanding what happens when 
we assume that f(n) = g(ni,n 2 , ■ ■ ■ ,njc), i.e. / depends on a finite number K 
of slow varying lattice variables nj such that g(rij ± k) = f(n ± kN 2 i) J2J). As 
we are mainly interested in applying in Section [31 these results to the lpKdV 
CJ, we need to know what happens to the function g when the function / is in 
the point n + 1. One needs to get explicit expressions for f(n + 1) in terms of 
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g(ni, ri2, ■ ■ ■ , tik) on different points in the ri\, n^, . . . ,rik lattices. At first let 
us study the case, considered in [12] when we have only two different lattices, 
i.e. K — 1. Using the results obtained in this case we will then consider the 
case corresponding to K — 2. The generic case will than be obvious. 
In the case of one variable we can use the result contained in [10]: 

00 i-l 

A k H g(n 1 ) = Y / -jP(hk)Ay(n) (3) 

i—k 

where H is any one of the possible Hj introduced before and Ajjg(nx) — 

12i=o(~ l) fc ~' (i).9( n i the & -variation formula obtained using a two-points 
forward difference scheme. The coefficients P(i, k) are given by 

with Sj , &l Stirling numbers of the first and second kind respectively. A table 
with the coefficients P(i,k) for (i, k) < (6,6) is contained in ref [10]. 

Eq. 10) allow us to express a difference of order k in the lattice of spacing 
H in terms of an infinite number of differences on the lattice of spacing h. To 
get an approximate solution we have to truncate the expansion at the r.h.s. of 
cq. by requiring a slow varying condition for the function f(n). We will say 
that the function f(n) is a slow varying function of order p if A^ +1 /(n) = 0. 
A slow varying function of order p is a polynomial of degree p in n [5] . For a 
function of order p — 1 , eq. J3J reduces to 

A H <?(ni) = N 2 A h f(n). (4) 

Dividing eq. (0} by h and taking the limit as h — > 0, with x = hn and x\ — 
n\Hi = mhN finite, we get d ^j^ = e dg }^ - ■ In the case p = 2, we get: 



A 2 H g( ni ) = N 4 A 2 h f(n), (5) 

V2 
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N 2 ( N 2 — 1 "1 

A H g( ni ) = N 2 A h f(n) + ^- ^/(n). (6) 



For p = 3 we have 

Aftf(m) - iV 6 A|/(n), 
A| ff (ni) = iV 4 A2/(n) + iV 4 (7V 2 - l)A 3 ,/(n), 
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and for p = 4 we have: 

AUim) = N s Afj(n), (7) 
A%g( ni ) = N^lf{n) + ZN *^ - l) A*/(n), 
A^(n x ) = iV 4 A 2 /(n) + 7V 4 (iV 2 - l)A£/(n) + 
+[ljV 6 (iV 2 _ i) _ liiv 4 (7V 2 - l)]A 4 /(n), 

A„g(m) = N 2 A h f(n) + ^N 2 (N 2 - l)A 2 /(n) + 

+ ±N 2 (N 2 1)(N 2 2)Alf(n) + ±N 2 (N 2 1)(N 2 2){N 2 - 3)A 4 /(n). 
From if /(n) is a slow function of order 1, f(n+ 1) reads: 

/(n + 1) = 5(n x ) + ^[ff(ni + 1) - <?(m)] + o(-^j), (8) 
while, if the function f(n) is a slow varying function of order 2, /(n + 1) is given 

by 

/(n + l) - ff ( ni ) + _^[_ s ( ni + 2)+4 5 (ni + l)-35(ni)]+ (9) 

+ i9(ni + 2) - 2 5 (m + 1) + g(m)] + o(^). 

When the function /(n) is a slow varying function of order 3, /(n + l) is given 
by 

/(n+l) = ff(ni) + ^[2 5 (ni+3)-9ff(ni + 2) + 13fl(ni + l)-6fl(ni)] + 
+ ^ i [-9(n 1 +3)+4 9 (n 1 +2)-5 9 (n 1 + l)+2 9 (n 1 )]+ (10) 

+ g^e + 3 ) - 3 5(«i + 2 ) + 3 s(™i + 1) - + K^g). 

and, when the function f(n) is a slow varying function of order 4, is given 

f(n + l)=g( ni )+ (11) 
+T2lvs[-30(n 1 + 4) + 16ff(m + 3) - 36<?(ni + 2) + 48 5 (n! + 1) - 25 S (m)] + 
+ 2ik[ n ff("i + 4) - 56 5 (nj + 3) + 114 5 (m + 2) - 104 S (m + 1) + 35 9 (n x )] + 
+ T2^[-35(«i + 4) + 14 5 (m + 3) - 24 S (m + 2) + 18 5 (m + 1) - 5g(ni)]+ 
+ 24^[5(«i + 4) - 4 fl (m + 3) + 6.g(nx + 2) - ig{n x + 1) + fl (m)] + o(^nr). 

In Section [21 we consider the reduction of an integrable discrete equation 
and will be interested in obtaining from it integrable discrete equations. It 
is known [20] that a scalar differential difference equation can possess higher 
conservation laws and thus be integrable only if it depends symmetrically on 
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the discrete variable, i.e. if the discrete equation is invariant with respect to the 
inversion of n. So we will choose asymptotic discrete formulas which contain 
both /(n±l). The results contained in eq. J3J| do not provide us with centralized 
formulas. To get symmetric formulas we need to take into account the following 
observations: 

1. Eq. © is valid also if H and h are both negative; 

2. For a slow varying function of order p, A?/(n) = A^/(n + I) for any 
integer number I. 

Using these observations, from eq. JHJ we get 

f(n - 1) = g{n x ) + ^[g(ni - 1) - <?(m)] + o(j^), (12) 
and in place of eq. © we have 

/(n + 1) = 9(m) + ^[g(n 1 + l)-g(n 1 -l)]+ (13) 
1 -\g{n x + 1) - 2 fl (m) + g(m - 1)] + o(-^), 



2iV 4 TV 6 ' 

when the function f(n) is a slow varying function of order 2. To get eq. I|13[) 
we have to write, using Observation 2, eqs. (0EJ) in the form 

A 2 H g( ni ) = iV 4 A 2 /(n), (14) 

n 2 (n 2 - ^) - 

A ffff (n x ) = JV 2 A fc /(n) + — ^ W/(*)> (15) 

where by the symbol A we mean the centralized version of the difference oper- 
ator. In the case of eq. ifUlfT^ A 2 /(n) = f(n + 1) - 2f(n) + f(n - 1). 

When f(n) is a slow varying function of odd order we are not able to con- 
struct completely symmetric derivatives using just two points forward difference 
formulas and thus /(n + 1) and f{n — 1) will never be expressed in a symmetric 
form (see for example eqs. (0 I12(l 'l. In the case when the function f(n) is a 
slow varying function of order 4, we have to rewrite formulas Q using both 
observations introduced before. We have: 

A^(m)=iV 8 A£/(n), (16) 

A^(nx) = N*Alf(n) + = 1} A*/(n), 

A 2 H g(m) = iV 4 A 2 /(n) + N\N 2 - l)A|/(n) + 

+ [^ 6 (iV 2 - 1) - ^N\N 2 - l)}AU(n), 

A„g(m) = N 2 A,J(n) + ±N 2 (N 2 - l)(N 2 - 2)A|/(n) + 

+ IjV 2 (JV 2 _ i )A 2 /(n) + i.iV 2 (7V 2 - l)(iV 2 - 2)(7V 2 - 3)Aj[/(n). 
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where, using Observation 2, we have written A|/(n) as h\f(n) — 2f(n + 2) — 
7f(n + 1) + 9f(n) — hf(n — 1) + f(n — 2). Inverting formulas l|16fl we have: 
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f(n + 1) = g(m) - ^ [g(m + 2) - 8 5 (m + 1) + 8 5 (m - 1) - 5(711 - 2)] (17) 
~24^4 + 2 ) - 16 5(«i + 1) + 3 °.9(«i) - 16.g(«i - 1) + S (ni - 2)] 
+ TTtW[5("i + 2 ) - 2 5("i + 1) + 2.g(m - 1) - ff(ni - 2)] + 



Let us pass now to the case of functions of multiple variables, i.e. when 
g = g(n%,n2, ■ ■ ■ , tik)- If n%, ri2, ■ ■ ■ , njt were completely independent discrete 
variables than any operation on one of the variables will not reflect on the other. 
If, however we have f(n) — 5(711, 712, ■ ■ ■ , Uk ) and rij = j^j, then any shift of 
n will reflect on all variables m, 712, ■ ■ ■ , n-^. We will consider later the case 
when the multiple variables are independent, i.e. the partial difference case. 

Let us consider here in all details the case of K — 2, which is the case we will 
need in the Section^ f(n) = 3(^1,^2) and we are looking for a representation 
of f(n + 1) in terms of <?(ni, 712) and its shifted values. The resulting formulas 
will depend in a crucial way on the slow varying order of the function f(n) with 
respect to n\ and 112- If the variation of both variables has to appear in the 
expansion of f{n + 1) than the slow varying order with respect to n\ must be 
greater than that of ri2- f(n) cannot be a slow varying function of order 1 in 
ni as in this case f(n + 1) will have no variation in n<i. If fin) is a slow varying 
function of order 2 in n\ than it can be either of order 1 or 2 in n%. In both 
cases the obtained formula will be valid up to order jh, but, in the first case the 
obtained expression will not be symmetric in ri2- When f(n) is a slow varying 
function of order 2 in n\ and of order 1 in 112, taking into account eqs. 1)151 Bjl 
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[g(ni + 2) - 4g(m + 1) + 65(711) - Ag(m - 1) + 5(711 
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and the observation given before, we have: 

fl(ni + l,n 2 ) = g(n 1: n 2 )+N 2 [f{n+l,n 2 ) - f{n,n 2 )}+ (18a) 
N 2 (N 2 - 1) 



-[f(n + l >na ) - 2/(n, n 2 ) + f(n - 1, n 2 )], 
= g(m, n 2 ) - N 2 [f{n + l,n 2 ) - f(n, n 2 )}-\ 

'-[fin + l,n 2 ) - 2/(n, n 2 ) + /(n - 1, n 2 )], 



2 

ff(m - 1, n 2 ) = ff(ni, n 2 ) - iV 2 [/(n + l,n 2 ) - /(n, n 2 )]+ (18b) 
iV 2 (7V 2 - 1) 



2 

g(n!,n 2 + 1) = ,g(ni, n 2 ) + iV 4 [/(ni, n + 1) - /(rii, n)], (18c) 

ff(ni + l,n 2 + 1) = g(ni, n 2 + 1) + iV 2 [/(n + 1, n 2 + 1) - /(m,n 2 + 1)]+ 

(18d) 

iV 2 (A^ 2 - 1) 

+ — ~[f{n + 1, n 2 + 1) - 2/(m, n 2 + 1) + /(n - l,n 2 + 1)], 

flf(m - l,n 2 + 1) = g(m, n 2 + 1) - iV 2 [/(n + 1, n 2 + 1) - /(m,n 2 + 1)] + 

(18e) 

N 2 (N 2 - 1) 

+ y — -[/(n + 1, n 2 + 1) - 2/(m, n 2 + 1) + f(n -l,n 2 + 1)], 

where we took into account that on the unshifted point n = ri\ = n 2 , f(n) = 
f(n,n 2 ) = f(ni,n) — f(n,n) — g(ni,n 2 ) and that f(n+ l,n 2 + 1), which is 
appearing in eqs. I|18dlll8cl) . is obtained from g(ni,n 2 + 1), given by eq. I|18c|) . 
by substituting n\ by n + 1. In such a way the 5 variables on the l.h.s. of eqs. 
l(T%|> are expressed in terms of the 6 variables f(n), f(n+ 1) = f(n + 1, n + I), 
/ (ni, n+1) = /(ra, n + 1), /(nil, n 2 ) = /(n±l,n), f(n— 1, n+1) appearing on 
the r.h.s. of eqs. (|18fl . Let us notice that to get a coherent number of equations 
with respect to the unknowns we had to consider also eqs. I|18bl I18e|) which 
involve n\ — 1. We can invert the system Ijl8(l and get: 

f(n + 1) = g(n u n 2 ) + wrnbini + l,n 2 ) - g{ni - l,n 2 )] + (19) 

+-^4 [g{ni, n 2 + 1) -g(ni,n 2 )] + 

+ - !> "2) - 2g(n 1 ,n 2 ) + g(m + 1, n 2 )]} + o(-^ ). 

In the continuous limit, when n — > oo and /i — > in such a way that nh — > a;, 
cq. <|19[1 will give 

/,x(a) = tg >Xx (x x ,x 2 ) + e 2 g yX2 (x 1 ,x 2 ). (20) 
When f(n) is a slow varying function of second order in both variables, in 
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place of eq. (|T%)l we have: 

g( ni + l,n 2 ) = g( ni ,n 2 ) + N 2 [f(n + l,n 2 ) - f(n,n 2 )}+ (21a) 
N 2 (N 2 - 1) 



g(m - l,na) = ff(ni, n 2 ) - JV 2 [/(n + l,n 2 ) - /(n, n 2 )]+ (21b) 

iV 2 (7V 2 - 1) 

L [f(n + l,na) - 2/(n, n a ) + /(n - 1, na)], 

5(n 1 ,7i 2 + l)= 5 (rii ! 7i 2 )+iV 4 [/(n 1 ,n+l)-/(n 1 ,n)]+ (21c) 
iV 4 (7V 4 - 1) 

+— ^ L [f(ni,n + 1) - 2/(m,n) + /(m,n - 1)], 

$(ni,na - 1) = $(ni,na) - N 4 [f(m,n+ 1) - /(m,n)]+ (21d) 
7V 4 (iV 4 - 1) 

+ -[/K,n + 1) - 2/(m,n) + /(m,n - 1)], 

<?(m + 1, n 2 + 1) = n 2 + 1) + iV 2 [/(n + 1, n 2 + 1) - /(m, n 2 + 1)] + 

(21c) 

[/(n + 1, na + 1) - 2/(m, n 2 + 1) + /(n - 1, n 2 + 1)], 

= ff(ni, n 2 + 1) - iV 2 [/(n + 1, n 2 + 1) - /(m,n 2 + 1)] + 

(21f) 

+ JV V " — \f(n + 1, n 2 + 1) - 2/(n l5 n 2 + 1) + /(n - 1, n 2 + 1)], 

= g{n l7 n 2 - 1) + AT 2 [/(n + 1, n 2 - 1) - /(m,n 2 - 1)] + 

(21g) 

[f(n + 1, n 2 - 1) - 2/(n l5 n 2 - 1) + /(n - 1, n 2 - 1)], 

= 5 (m,n 2 - 1) - 7V 2 [/(n + l,n 2 - 1) - /(m,n 2 - 1)] + 

(21h) 

+ - VJ< 2 — [f(n + 1, na - 1) - 2/(m, n 2 - 1) + /(n - 1, n 2 - 1)]. 

In this case the 8 variables on the l.h.s. of eqs. (|21|l are expressed in terms of the 
9 variables /(n) — f(n,n), /(nil) = /(nil, nil), /(ni,nil) = /(n,nil), 
f(n i 1, n 2 ) = f(n i 1, n), f(n — 1, n + 1) and /(n + l,n — 1) appearing on the 
r.h.s. of eqs. (|21|l . One can invert the system l|21|l and gets: 

/(n+1) = g(n 1 ,n 2 ) + ^\ 9 (n 1 + l,n 2 )-g(n 1 -l,n2)]+ (22) 

[fli(ni - 1, n 2 ) - 2g(m,n 2 ) + g{n\ + 1, n 2 )] + 



iV 2 (7V 2 


-1) 


2 




- l,n 2 


+ 1) 


iV 2 (iV 2 


-1) 


2 




+ l,n 2 


-1) 


N 2 (N 2 


-1) 


2 




- l,n 2 


-1) 


7V 2 (iV 2 


-1) 



2iV 4 

+ ^4 [g(ni,n 2 + 1) - g(ni,n 2 - 1)] + o(j^). 

In the continuous limit, when n — ► oo and h — ► in such a way that n/i — * .t, 
eq. (J52J will give eq. f2"U|) . 
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Let us write here just the final result when f(n) = g(ni, n 3 ) is of order 4 in 
the variable ri\ and of order 2 in 77,3. We have: 

/(n + l)= 5 (n 1 ,n 3 )- (23) 
- \9i n i + 2 , n 3) - 85(711 + 1, n 3 ) + 8#(m - 1, n 3 ) - g-(m - 2, n 3 )] + 

+ ^^4 [5(«i + 2,n 3 ) - 16g(ni + l,n 3 ) + 30g(ni,n 3 ) - 16g(m - l,n 3 ) + 

+ 5 (m - 2, n 3 )] + ^6 (m, n 3 + 1) - g(ni,n 3 - 1)] - [g{ni + 2, n 3 ) 

-2g(m + 1, n 3 ) + 25(711 - 1, n 3 ) - g(ni - 2, n 3 )] + o(j^). 

In the continuous limit, when when n — > 00 and ft, — > in such a way that 
n/i — > x, eq. will give 

f,x(x) = £g, Xx (xi,x 3 ) + e 3 g yX3 (x!,x 3 ). (24) 

One can introduce constant parameters in the definition of n\, n 2 or n 3 in 
terms of n. For example we can write n\ = ^Mr and n 2 = n^ar- Mi and M; 
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cannot be completely arbitrary as ni, n 2 , n and iV are integers. In such a case 
eq. fH^I reads: 

Mi 

+ = g{ni,n2) + -^p[g(ni + l,n 2 )-g(ni- l,n 2 )] + (25) 

+ ^ls(™i,n 2 + l)-p(ni,n 2 )] + 

M 2 1 
+ ^?{[.9(^i - i)^) - 2g(ni,m)+g(ni + 1, ra 2 )]} + o( — ). 

In Section 13 we will apply these results to a partial difference equation. For 
the sake of simplicity from now on we write the independent variables as indices. 
For completeness, in the following we present the formulas for two independent 
lattices, n and to, and a function defined on them f n>m . As the two lattices 
are independent the formulas presented above apply independently on each of 
the lattice variables. So, for example, f n +i. m when the function / is a slowly 
varying function of order 2 of a lattice variable ri\ (see eq. <|13[1 ) will read: 

fn+l,m — <7m,m + 2]\[2 \8ni+l,m ~ g-ni-l.m } + (26) 

— 2g„ lim + g ni -i,m] + °(t7^)> 



and similarly for a variation with respect to to alone or to the case when we 
will introduce multiple lattices associated to n or to, when formulas 1)191 1221 l2*3*)l 
are to be taken into account. A slightly less obvious situation appears when we 
consider / n +i,m+i, as new terms will appear. We consider here just the case we 
will need later when n\ = -^r 1 , m\ — M £™ and to 2 = jpc- If / is a slow varying 
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function of first order in 7712 and of second order in both n\ and mi, /n+i.m+i 
reads: 

Mi 

] + (27) 

T 27V 2 ^" l!?Ttl + 1 ' m2 ~ 9n 1 ,m 1 — l,m 2 \ T 

"T 2iV 4 ^ ni + 1,mi '" 12 ^~ 9ni — l,m 1 ,m 2 ^9n 1 ,m 1 ,m 2 \ T 

M I r 

27V 4 l-f™l,™l + l>™2 9ni,mi — l,m.2 < ^9ni ,m\ ,m 2 \ 4" 
MiM 2 r 

4./V 4 l5' r *i+l ) 'Wi+l,ma "r — l,mi — l,m 2 ffni+l!»r*i — l,mj 
— AW — l.mi+ljHis] + ^4 [Stii.mi ,m 2 + l — 9n 1 ,m 1 ,m 2 ] + °^7V^^' 

As one can see in its fifth and sixth lines, eq. (|27|l contains extra terms involving 
variations in both the mi and the ni lattices. 



3 Reduction of the lattice potential KdV 

In this Section we apply the results presented in Section 2 to the case of the 
lattice potential KdV QJ. By expanding the left hand side of eq. JU, one 
separates the linear and nonlinear parts: 

(j) - q)(u n+ i <m+1 - U n , m ) + {? + q){Un+l,m ~ "n,ro+l) = (28) 

l)(Un+l ,m+l ^n^m)' 

This equation involves just four points which lay on two orthogonal infinite 
lattices and are the vertices of an elementary square. 
Let us solve the linear equation 

F = (P - ?)("n+l,m+l - u n,m) + (p + <?)( u n+l,m ~ «n,m+l) = °- (29) 

The discrete Fourier transform [5] will reduce the solution of the Partial Differ- 
ence Equation (PAE) i|29|l to that of an ordinary difference equation. Defining: 

u n>m = tt~ f v m (z)z n ~ 1 dz, (30) 
2tti J Cl 

v m (z) = Un < rn ( 31 ) 

n=— 00 

where Ci is the unit circle, we reduce equation l|29|l to the following first order 
equation for v m (z) 

v m +i(z)[(p - q)z - (p + q)} - v m (z)[(p - q) - z(p + q)} = 0, (32) 
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whose solution is given by 

v m (z) = 



(p-q)- z(p + q) 



v (z). 



(p-q)z -(p + q)- 

Given any initial condition u n fl the general solution of 129J1 is given by 



2m 



i E 



«i,0 



] = -CO 



Ci 



(p-q)- z(p + q) 



(p - g)z -(p + q) 



z n -i- l dz. 



(33) 



(34) 



Eqs. (|30l I34|l can be rewritten in a more natural way (from the continuous point 
of view) by defining 



n = 



(p-q) - z(p + q) 



(p - q)z - (p + q) 



(35) 



In such a case eq. (|3()t is just the standard Fourier transform and the solution 
<|34|) is just written as a superposition of linear waves. The dispersion relation 
for these linear waves is given by lo — u)(k) and reads: 



cj = —2 arctan 



— tan ( — 
2 



(36) 



In the following, however, to avoid too complicate formulas we express the 
solutions of the linear equation in term of z and f2. 

The lpKdV is an integrable equation of the same category of the KdV [3] as 
it possesses a Lax pair [15] which can be obtained by requiring that the model be 
consistent around a cube. So, as from KdV we get by multiple-scale reduction 
the NLS [21], the same we may expect here [2]. To get an integrable discrete 
equation we expect a resulting discrete equation which is somehow symmetric. 
At least when ht — ► with t = mh t the differential difference equation we 
obtain must be symmetric in terms of the inversion of rij [20], i.e. if it depends 
on rij + k it will depend in the same way on nj—k- So in the transformation of the 
discrete dependent variables we prefere to use formulas i|13l H7|) for the space 
variables, while considering the lowest possible approximation for the discrete 
time variable ©. So, in the multiple-scale expansion, as we do not need to 
have the discrete time variable appearing in a symmetric way, we use eqs. (I19J) . 

Taking into account eq. 134|) , we consider a wave solution of eq. ( I29fl given 

by 



E n 



(37) 

Eq. H37|) solves F — if u> is given by eq. I|3()|) . Then we look for solutions of 
eq. H28fl written as a combination of modulated waves: 



+ oo 



^ t e S "4 , n 1 ,m 1 ,m- 2 (En,m) + 
s=0 

+ OG 



(38) 



s=0 
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where ftlmi.ma are slowly varying functions on the lattice and e = N~ 2 . 
By a we mean the complex conjugate of a so that, for example, E n ^ m = 
e -i(kn-u(k)m) _ (E^m) -1 , and the positive numbers (3 S are such that /3o = 1 
and Pj = j. The discrete slow varying variables n\, mi and m 2 are defined in 
terms of n and m by 



A r2 



Mi' 

TV 2 
M? 
TO 2 iV 4 . 



mi- 



(39) 



Eq. 1)39(1 is meaningful if Mi and M2 are divisors of N 2 . 

Introducing the expansion l(38|l into eq. 1(28(1 and picking out the coefficients 
of the various harmonics (E n ^ m ) s we get a set of determining equations. For 
s = l, having defined tp^ = ip, we get at lowest order in e: 



(q-p)(l-zn)-(p + q)(n-z) 



0. 



(40) 



which is identically solved by the dispersion relation 136(1 . At e 2 we get a linear 
equation 



M 2 
-Mi 



(p - q)zVL -{p + q)Sl) 



V^ni ,mi +1,7712 V^Tll ,7711 — 1,7712 

—1,7711 i m 2 



= 



(p - q)zQ + (p + q)z V'm+i.mi.m 
whose solution is given by choosing 

^ni,mi,m2 — 4 > n2,rti2 

where 

712 = n\ — m\. 

The solution 1(42(1 is obtained by choosing the integers Mi and M 2 as 
Mi = S fi[(p-g)«-(p + g)], 
M 2 = S z\(p-q)n + (p + qj\, 



(41) 

(42) 
(43) 
(44) 



where S is an arbitrary complex constant. Let us notice that also n 2 = n± +m± 
solves eq. lj4*3*)l by an appropriate choice of Mi and M 2 . Moreover, = u t k, 
the group velocity. As Mi and M 2 are integers, not all values of k are admissible 
as uj^k must be a rational number. 



At e we get a nonlinear equation for <p n2 m2 
A Q ) _<a(°) 

T / n 2 + l,rn2 Tn 2 -l,m2' 

) + C 2 ((f>n 2 +2,m 2 + 4>n 2 -2,m 



which depends on ip, 



(2) 

ri2 ,m 2 



and 



2<^n 2 ,m 2 ) "H 



(«) 



n 2 ,m 2 Vrn 2 + l,m 2 



(45) 
<// 0) 1 + 



+C5V'r? 2 ) ,m 2 ^ 



m 2 Vn 2 ,m 2 



= 0, 
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where 
c\ 

(-2 
C3 
C4 



[(p-q) -z(p + q)], 



S 2 z 2 pq(p - q) 



-25 zpq(p-q) 



(p + q)z-(p-g) -} 2 

(p + q) - z(p-q)\ 
Av + q) - (P - q)] [{z 2 + !)(? + «)- 2(p - g)] 



c 5 = -2 



25V-g 2 )(z 2 -l) 



[(P + 9) - z(p-q)} 2 

Z+l l 2 



(p + q) - ~Av - q) 

(z + l)(z 2 -l) ^ 



z[(p + q)z - (p - q)} l(p + q)- z{p - q) 



The lowest order equations for the harmonics s = and s = 2 appear at e 2 and 
give: 



"«8+l,in! rn 2 — l,m2 



2|0„ 



(1 + z) 2 



5pz[(p + g)z - (p-?)]' 



2 1 + z 

u n 2 ,m 2 } c, ft \ 

2p(l - z) 



Taking these results into account the nonlinear equation (|45|l for 



^n 2 ,m 2 + l i vn 2 ^ra 2 



C\(k) 



f ) n 2 +2,m 2 + 4>n 2 -2,; 



20r, 



(46) 
(47) 

reads: 
(48) 



-C 2 (k) 



^?12 + I,m2 T*r|.2 — 1,171-2 

-20 n2 , m2 +C 3 (fc) 

n 2 ,m 2 \ ( f > n2 ,ni2\ ? 

where Cs(fc) is a real coefficient given by 

C 3 ( fe )=2_ sin( fc )(l + cos( fc)) 2 g (p 2 ^ 2 ) 

p f(p 2 - 9 2 ) 2 (cos (fc)) 2 - 2 (p 4 - g 4 ) cos (Jfe) + (p 2 + q 2 ) 2 ) 

The coefficients Ci(fe) and Ci{k) are complex and depend on 5. They are: 

(z(p + q)-(p- q)) 



C\(k) = iqpz 2 S 2 {p-q) 



{{p-q)z- (p + q)Y 



C 2 {k) = -2iqpzS 2 (p-q) 



((p + q)(l + z 2 )-2(p-q) z) 
{{p-q)z- {p + q)) 2 



(50) 



(51) 



Eq. I|48|l is a completely discrete and local NLS equation depending on the first 
and second neighboring lattice points. At difference from the Ablowitz and 
Ladik [1] discrete NLS, the nonlinear term is completely local. 
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4 Discussion of the results and conclusive re- 
marks 



The choice of the order of slowlyness is essential in defining the points involved in 
the resulting equation. In our calculation of the multiple-scale reduction of the 
lpKdV we choose to use the minimum number of points in the various lattices 
introduced. Here we started from just four points and got a scheme which 
involves six points. Moreover while the starting initial problem is defined on a 
staircase, eq. (|48|) is defined on a line. By choosing slow varying functions of 
higher order, essential for example to get the higher order terms in the expansion 
necessary to go beyond the NLS [4] even at the lowest order we would get a 
nonlinear difference equation involving many more lattice points. This seems to 
be a peculiarity of the multiple-scale expansion on the lattice. 

This work open a research field of great interest both for the possible math- 
ematical results and for the physical applications. To show this one presents in 
the following a detailed list of open problems and remarks on which work is in 
progress: 

1. Prove the integrability of eq. (|48|l by reducing the Lax pair of the lpKdV 
or by constructing its generalized symmetries. 

2. Eq. (|48|) is invariant with respect to time translation. One can thus reduce 
it with respect to this Lie point symmetry and get: 



!J n 2 +2 



u n 2 -2 



2</>„ 2 + di 0„ 2+ i + 0„ 2 _i - 20„ 



(52) 



12 <Pn 2 \Pn 2 \ 



= 0. 



One would like to show that equation (|52|) possesses the Painleve prop- 
erty. On the lattice this is given by the singularity confinement [7, 16] or 
algebraic entropy [9]. 

Eq. I|48|) has a natural semi-continuous limit when m2 — > oo as H2 — > 
in such a way that t<i — m,2H2 is finite. In such a case eq. ()48f) reduces to 
the nonlinear differential difference equation 



dU 



ei 



^2+2 + 07l 2 -2 - 20„ 2 

K2+1 + 4>n 2 -i ~ 2</>n 2 



+ 



(53) 



?3 fi^lfin 



If eq. H48[) is integrable than eq. H53[) should be an integrable multiple- 
scale reduction for differential difference equations [11-13,17] like the Toda 
lattice. 

Do the multiple-scale reduction of other integrable lattice equations like 
the time discrete Toda lattice, the lattice mKdV or the discrete sine- 
Gordon equation and see what one gets. One could obtain eq. (|48|1 but, 
maybe some other integrable lattice NLS like equation can be obtained. 
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5. Do the multiple-scale reduction of the discrete Burgers equation 

On,m+l — 1 7 7 [7—7 ' (0 4 J 

J- ~r On,m°n — l,m ~r a "n-l,m 

and get a discrete Eckhaus equation [2]. 

6. Apply the reduction technique to some nonintegrable equation of physical 
interest like for example those obtained in the case of discrete phenomena 
in liquid crystals [6] and obtain approximate theoretical solutions of the 
physical result. 
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